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SUMMARY 

Combustion instability phenomena are examined theoretically 
using a small amplitude, linear analysis. Conservation equations 
are solved for finite axial gas velocities in the combustion chamber. 
Transverse waves with axial variations were considered. Concen- 
trated combustion near the injector face was assumed. The gen- 
eral wave equations are solved for specific boundary conditions at 
the end of the chamber (nozzle entrance), walls of the chamber 
(acoustic liners) and the injector (combustion zone). 

Calculated curves are presented to illustrate how Mach number, 
chamber length to radius ratio, nozzle flow response, and acoustic 
liners influence the combustion response required to obtain neutral 
stability or various growth rates. These curves demonstrate the 
importance of the perturbations in combustion rate that are 90° out 
of phase with the pressure perturbations. The variations in pressure 
amplitude and phase, or time displacement, with axial position are 
presented to illustrate how the combustion process can influence the 
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wave characteristics. For example, it is shown that a combustor 
can have oscillations over a range of frequencies that are dependent 
on the combustion process. 

Results of stability tests with a 20 000-pound thrust hydrogen- 
oxygen combustor are briefly analyzed with the aid of the calculated 
neutral stability limits. The comparison illustrates how the com- 
bustion response calculations can be used to predict stability or what 
acoustic liner is required to stabilize a combustor. 

INTRODUCTION 

Recent combustion instability experience has demonstrated several 
interesting phenomena that have not been described theoretically. 

Most interesting of these is that acoustic absorbing chamber walls 
change stability limits and frequencies of the instability (see ref. 1). 
Detailed pressure measurements have shown (refs. 2 and 7) that 
waves travel as transverse acoustic waves but they do not have con- 
stant amplitude and phase with length. Other data indicate that the 
frequencies of the instability do not always correspond to those calcu- 
lated for pure acoustic modes. These theoretically unexplained re- 
sults motivated this study. 

The variations in the transverse wave properties with length 
suggest that modifying the theoretical approach of ref. 3 to include 
axial variations would be successful. The data with acoustic liners 
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indicated that the boundary conditions at the walls would be controlled 
by the liner. Therefore various boundary conditions would have to 
be included in the theory. Recent publications (refs. 4 and 5) 
wherein the combustion process is characterized by a response func- 
tion (perturbation in burning produced by a perturbation in pressure) 
suggest that a term like a response function be used to generally 
characterize the combustion process. Since most combustors operate 
at high gas velocities or Mach numbers within the chamber it also 
seemed necessary to allow for high Mach numbers. 

All of these aspects were incorporated in a theory to obtain a 
relatively simple equation to characterize the stability characteristics 
of a system. Calculations were then performed to illustrate how 
various design parameters influence the stability characteristics. 
Experimental stability characteristics of a 20 000-pound thrust hydro- 
gen oxygen combustor are compared to the calculated results to il- 
lustrate how to use the theoretical results and the validity of the theory. 

THEORY 

For this analysis it is assumed that combustion is concentrated in 
a zone immediately adjacent to the injector face. It is also assumed 
that during unsteady operation, small amplitude, three-dimensional, 
irrotational waves propagate through the chamber. These waves must 
satisfy the boundary conditions of the chamber. The waves produce 
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a perturbation in the combustion rate at the injector face which is 
determined by the impedance of the combustion process. The waves 
also produce a perturbation in flow at the end of the cylindrical cham- 
ber, or entrance to the convergent nozzle, which is determined by 
the impedance of the nozzle. Similarly, if the walls of the cylindri- 
cal chamber have acoustic absorbing characteristics, the waves will 
produce a perturbation in velocity at the wall which is also determined 
by an impedance. 

The equations of motion of an inviscid compressible fluid, which 
describe the waves in the chamber can be written in nondimensional 
form {p, P and T are nondimensionalized by their average values, 



yields (see ref. 3) 
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Now $ and P can be expanded in a series in powers of some ampli- 
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Since V5>(0) = M then eq. (2) yields the following terms containing 
e to the first order. 
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If it is assumed that the solution of equation^ is periodic in time and 
separable in the Z-, r-, and 9 -coordinates the following solutions 
are obtainable. 
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In eqs. (5) to (10), the terms n, m, co, Bp B 2 , and C carry 
the following connotation. 

n number of waves that exist in the 6 direction and must 

be an integer 

co complex number, the real part specifying the frequency 

of the oscillation and the imaginary portion specifies the 
growth or decay rate with time 
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m complex eigenvalue of the radial differential equation as 

determined by the wall boundary conditions 
Bf, B 2 complex numbers whose real values specify the speed at 

which a wave travels axially up or down the chamber. The 
imaginary part specifies the decay in amplitude with 
length 

C complex number which characterizes the waves at the noz- 

zle end of the chamber 

A flow response is assumed to characterize the boundary conditions 
at the nozzle and injector and a velocity response at the wall as given 
by: 
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Substituting eqs. (7) and (8) in eq. (11) to determine the nozzle response 
at the nozzle end of the chamber (Z = 0 ) and solving for C: 
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Substituting eqs. (7) and (8) in eq. (11) to determine the combustion 
response N at the injector end of the chamber (Z = -L): 
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Substituting eqs. (7) and (9) in eq. (12) to determine the wall response 
K at the chamber wall (r = 1): 
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The amplitude (P max ) of the wave, determined from the product of 
the pressure and its complex conjugate, and the angular position (/ 3 ) 
of maximum amplitude at any time, and radial or axial position is 


given by: 
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Equations (6), (13), (14), and (15) were solved on an IBM 7094 
computer for various values co, m, M, L/R, and G and for n = 1 
(first transverse mode) to obtain combustion responses (N) and wall 
responses (K). For hard walls or K « oi, a value of m of 1. 84118 
is used. For soft wall (K + a), since K varies with chamber length, 
a trial and error technique was used to determine which value of m 
produced the desired average value of K over the full chamber length. 
Pressure profiles were determined from equations (6), (13) and (16) to 
determine maximum amplitude and phase of the wave relative to the 
injector end of the chamber as a function of length. 

RESULTS AND DISCUSSION 

Stability calculations were performed to examine and illustrate 
the influence of various parameters (i. e. , decay coefficient, length 
to radius ratio, Mach number, nozzle response and wall response) on 
stability and wave characteristics. A combustor with L/R = 2.7, 
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Mach number = 0. 33, nozzle impedance of Re(G) = 0. 9166 and infinite 
wall impedance was used as a reference condition (to illustrate how 
varying one parameter at a time from the reference combustor changed 
the stability and wave characteristic). The reference combustor condi- 
tions were selected to match those of a 20 000-pound thrust hydrogen 
oxygen engine extensively tested at the Lewis Research Center (ref. 6). 

Stability lines for various growth or decay coefficients are shown 
in figure 1. The neutral stability line (waves will not grow or decay) 
shows that for a combustion response of Re(N) = 0 and Im(N) = -0. 7 
a transverse wave of 0. 9 the frequency of a pure acoustic wave exists. 
As the combustion response is increased the frequency increases and 
approaches the pure acoustic frequency when the combustion response 
is near the nozzle response; Re(G) = 0. 9166 and IM(G) = 0. 0. With 
the real part of the combustion response above the value for neutral 
stability the wave grows with time. A decrease in response results 
in a wave that decays. Decreasing the imaginary part of the combus- 
tion response from the neutral stability value results in a wave that 
grows with time. Decreasing the real part of the combustion response 
or increasing the imaginary part improves stability. To illustrate how 
various parameters influence the stability characteristics of a system, 
only the neutral stability line will be used. 
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Influence of Combustion Chamber Design on Neutral Stability 

The influence of Mach number on neutral stability and wave pro- 
files are shown in figs. 2 and 3. With low Mach numbers the stability 
line is almost vertical. The stability limit is therefore uniquely de- 
fined by the real portion of the combustion response. Frequencies 
with low Mach numbers are very close to those of a "pure" acoustic 
mode in a chamber with hard walls and ends. At very low Mach num- 
bers the wave amplitude is constant with length and there is no phase 
displacement in the axial dimension. 

At higher Mach numbers the neutral stability lines are generally 
rotated about the value of combustion response equal to the nozzle 
response. This rotation is accompanied by an increase in the varia- 
tion of wave amplitude with length and phase shifting as shown in 
fig. 3. At low frequencies the amplitude at the nozzle end of the cham- 
ber leads the wave at the injector end (positive /3 angle). The ampli- 
tude variation lowers the real part of the combustion response needed 
to drive the wave. 

At higher frequencies the opposite effects are observed. Amplitude 
is highest at the nozzle and the wave at the injector end leads the noz- 
zle end. Generally it is found that when the real part of the combustion 
response is less than the real part of the nozzle response, the ampli- 
tude at the nozzle is less than at the injector. If the imaginary part of 
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the combustion response is less than the nozzle response the nozzle 
end of the wave leads the injector end. 

Increasing the Mach number also increases the range of frequencies 
over which the system can oscillate. Within the limits of combustion 
response shown in fig. 2 a Mach number of 0.01 has a frequency varia- 
tion of 0. 1%, while a Mach number of 0. 6 has a frequency variation of 
approximately 20%. Higher Mach numbers are also accompanied by 
several loops in the stability map as illustrated by the two loops shown 
for Mach number 0. 6. These loops correspond to the first transverse 
mode combined with the various longitudinal modes. Each time the 
curve crosses itself it corresponds to going to one higher longitudinal 
mode in the combined mode. For a Mach number of 0. 6 and at a re- 
duced frequency co/u> 0 = 1.26 the curve has crossed itself twice (two 
loops) and at this frequency the second longitudinal mode combined 
with the first transverse should be observed. The second mode is 
observed in fig. 3 for Mach number 0. 6 and w/o> 0 = 1.26 as there 
are two distinct minimums in the amplitude curve, corresponding to 
two pressure nodes. All the curves presented herein will form loops 
at higher frequencies. 

The influence of chamber length to radius ratio (L/R) on the 
neutral stability lines and wave characteristics is shown in figs. 4 
and 5. Decreasing the L/R generally moves the neutral stability line 
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down and to the right, thereby increasing the stable region of the sys- 
tem. Increasing the L/R produces more loops in the region where 
the combustion response is close to the nozzle impedance. Low L/R 
systems do not exhibit a significant phase shifting with length but do 
have appreciable amplitude variations as shown in fig. 5. The general 
characteristics of the influence of L/R are similar to those of Mach 
number as described previously. 

The influence of nozzle flow response on neutral stability lines is 
shown in fig. 6. Reducing the real part of the nozzle response moves 
the curves to the left and decreases the stable region. Increasing the 
imaginary part of the response eliminates the loop and generally moves 
the curve down and to the right, thereby increasing the stable region. 
Improved stability is achieved by increasing the real and imaginary 
parts of the nozzle response. 

Influence of Acoustic Liners 

The influence of an acoustic liner is shown in fig. 7. Various liners 
are shown with three having the same absorption coefficient, ot, which 
is generally used to characterize liners (ref. 1). For an absorption 
coefficient of 0. 1 the three liners gave about the same results, lower- 
ing the curves and shifting them to the right, increasing the stable 
regime relative to the hard wall system. The greatest improvement 
occurs at higher frequencies. Further increasing the absorption coef- 
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ficient to 0. 33 gave about three times the improvement achieved with 
a 0. 1 coefficient. 

Improvement in neutral stability lines with a 0. 1 absorption coef- 
ficient liner for different L/R chambers is shown in fig. 8. Sur- 
prisingly, the L/R does not significantly change the effect of the 
liner even though more surface to volume is available to absorb the 
wave with large L/R's. At large L/R the amplitude at the nozzle 
with a liner is much lower than without a liner. At small L/R the 
amplitude reduction is small. Thus, the increase in wall loss is 
compensated by a decrease in nozzle loss. 

Improvement in neutral stability lines with a 0. 1 absorption coef- 
ficient liner for different Mach numbers are shown in fig. 9. Mach 
number strongly controls the affect of an acoustic liner. At low Mach 
numbers the liner is extremely effective in increasing the stable 
regime. At higher Mach numbers the effectiveness of the liner is de- 
creased. The large increase in effectiveness at low Mach numbers is 
again due to the change in amplitude at the nozzle. For low Mach num- 
bers, the amplitude at the nozzle with a liner is greater than without a 
liner. Therefore the liner increased the loss through the nozzle as well 
as through the wall. At high Mach numbers, the liner decreases the 
nozzle end amplitude and losses. 
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Combining Combustion Proces s and System Responses 

The previous sections have shown what combustion response is re- 
quired at the injector end of a given combustor configuration to obtain 
various stability characteristics. The combustion process has a resp- 
onse which is determined by the injector, propellants, etc. A typi- 
cal response for vaporization of liquid oxygen droplets is described in 
ref. 4 and shown in fig. 10. The vaporization response is determined 
by the dimensionless time wr (co is frequency of oscillation and t 
is mean drop lifetime). At cur = 0 the real and imaginary parts of 
the vaporization response are zero. At infinite wt the vaporization 
response goes to Re(N) = -3. 93 and Im(N) =0.0. Also shown in 
fig. 10 is the combustion response required for neutral stability for 
the reference combustor configuration described previously . System 
stability characteristics for the vaporization process with the reference 
combustor are obtained when the combustion response line intersects 
the vaporization response line at the same frequency. The lines inter- 
sect when co t = 0. 2 and co/co^ = 0. 98 and when cot = 2.2 and 
co/co 0 = 0. 4. Therefore, neutral stability would be obtained for 
r = 0.2/co 0 and r = 5. 5/co 0 . Between these r values the system 
would be unstable as the vaporization response line would intersect 
combustion response lines of positive growth rates. Outside of these 
r values the system would be stable as the intersection would be with 
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lines of negative growth rates. Neutral stability is obtained with a 
combustion response of Re = 0. 5 and Im = -0. 03 and with Re = -2. 25 
and Im = -3.2. These real values are much less than the 0. 9166 values 
of the nozzle response. This illustrates why it is important to know the 
in-phase and out-of-phase portion of the combustion process response. 

Studies with the reference combustor have shown that stability was 
very dependent on the hydrogen temperature or density (ref. 6). To 
explain this effect a model was postulated (ref. 5) in which stability was 
influenced by flow response of the gaseous hydrogen. The response of 
the hydrogen flow and oxygen vaporization process was obtained by a 
flow weighted average of the response of the oxidizer and the hydrogen, 
or 

W w f 

N=w^N oxt ^N f (IT 

Curves of hydrogen response (Nj ) were calculated for different hydro- 
gen densities. Results of these calculations for the "Reference 1" 
system of ref. 5 are shown in fig. 11. Response of the Lq X is 
Re(N) = 0. 56 and Im(N) = -0. 02 obtained from fig. 11 at an oj j = 0. 2. 
The combustion response curve crosses the combined hydrogen flow 
and L qx vaporization response line at a hydrogen density of 0. 48 
which is within the 0. 4 to 0. 5 density limits obtained experimentally 
(fig. 5 of ref. 5). 
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Studies with the same combustion chamber showed (ref. 1, figs. 

29 and 30) that an acoustic liner with an absorption coefficient of 0. 08 
reduced the hydrogen temperature from 110° Rwith no liner to between 
95° and 60° R with the liner. This corresponds to a hydrogen density 
increase from 0. 48 to between 0. 59 and 0. 98. The combined hydrogen 
and oxygen response curve of fig. 11 intersect the 0. 1 absorption coef- 
ficient liner calculations at a density of 0. 90, which is in reasonable 
agreement with experiment. Better agreement would have been ob- 
tained if a larger cut had been assumed for the oxidizer. This also 
shows the importance of the imaginary part of the combustion response. 
Increasing the cut of the oxidizer mainly lowers the combined hydro- 
gen and oxygen response curve of fig. 11. Since without a liner the 
combustion response is almost vertical at the point of intersection a 
small change in cut of the Lqx would not influence the stability 
limits. With a liner added the influence would be very significant as 
the liner curve is not vertical. 

CONCLUSIONS 

A mathematical model to determine stability characteristics of 
rocket combustors with high Mach number flow and acoustic liners has 
been presented. Calculations performed with this model for the first 
transverse mode have shown many characteristics observed in actual 


combustors. 
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1. Variation in wave amplitude with axial position 

2. Phase or time displacement of the wave with length 

3. Strong influence of Mach number and length to radius ratio on 
stability limits 

4. Range of possible frequencies of oscillations with a given com- 
bustor configuration 

With this model it is possible to determine what acoustic liner is 
required to stabilize a combustor. Calculations have illustrated the 
importance of the perturbation in burning that occurs out of phase 
(90°) with the pressure oscillations. 

NOMENCLATURE 

A cross sectional area of combustor 

a speed of sound 

Bj, B 2 complex coefficients given by eq. (6) 

C complex coefficient given by eq. (13) 

D total derivative 

e exponential 

G complex number specifying nozzle response eq. (11) 

g gravitational constant 

Im( ) imaginary part of term in ( ) 

i unit complex = -yf-T 

J n Bessel function of Order n 
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K 

L 

M 

m 

N 

n 

P 

P 

max 

Q 


Re( ) 

r 

T 

t 

u 


W 

Z 


velocity response at chamber walls, eq. (12) or (15) 
length of chamber 
Mach number 

argument of bessel function 
combustion response eq. (11) or (14) 
number of pressure modes in 9 direction 
pressure 

amplitude of pressure oscillations 

constant of integration from integrating the momentum 

, y - 1 2 

equation = 1 + — ^ — M 

real part of ( ) 

radial dimension 

temperature 

time 

velocity vector 
velocity in tangential 
velocity in radial direction 
velocity in axial direction 
flow rate 
axial direction 



- 20 - 


a 


absorption coefficient = 


4Re(K) 


|Re(K) + l] 2 +[lm(k)] 


/3 angular position of maximum amplitude 

7 specific heat ratio 

e perturbation amplitude 

6 tangential direction 

p density 

t average vaporization time 

$ velocity potential V(<E>)=u 

co frequency of oscillation 

co Q frequency of oscillation = a x 0. 586/2R 

V del operator 

Subscripts: 


f fuel 

ox oxidizer 

t derivative with respect to time 


Superscripts: 

(1) perturbation quantity to order (1) of perturbation ampli- 

tude ( ) 


time average value 
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Figure 1. - Stability lines for different wave growth rates. Chamber length/ radius, 
2.7; Mach number flow, 0.33, Nozzle response, Re(G)- 0.9166, Im(G) - 0, hard 
wall chamber. 



Figure 2. - Influence of Mach number on neutral stability limits. Chamber 
length/radius, 2.7; nozzle response, Re(G), 0.9166, IM{G), 0, hard wall chamber. 
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(c) Chamber length/radius, 5.0 Figure 6. - Influence of nozzle flow response on neutral stability limits. Cham- 

ber length/radius, 2.7; Mach number flow, 0.33, hard wall chamber. 

Figure 5. - Pressure profiles with various length to radius ratio and frequencies. 

Mach number, 0.33; nozzle response, RE(G) * 0.9166, IM(G) - 0, hard walls. 
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Real part combustion response, Re(N) 


(b) Mach number, 0.33. (c) Mach number, 0.6. 

Figure 9. - Influence of an acoustic liner on stability limits with different 
Mach numbers. Length to radius ratio, 2.7. Nozzle response, RefG) 
-0.9166, Im(G) = 0. 



Real part combustion response, Re(N) 

Figure 10. - Combined response curves of oxygen vaporization 
with reference combustion chamber. Mach number, 0.33, 
nozzle coefficient, Re(G) * 0.9166, Im(G), hard wall. 
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Figure 11. - Combined combustion system response curves with hydrogen flow 
andLOX vaporization response. Mach number, 0.33; length to radius ratio, 
2.7; nozzle response, Re(G), 0.9166, Im(G) ■ 0. 
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